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Abstract 

We study the lattice N=l Wess-Zumino model in two dimensions and we construct a sequence 
p(G of exact lower bounds on its ground state energy density p, converging to p in the limit L — > oo. 
The bounds p^ L ' can be computed numerically on a finite lattice with L sites and can be exploited 
to discuss dynamical symmetry breaking. The transition point is determined and compared with 
recent results based on large-scale Green Function Monte Carlo simulations with good agreement. 
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I. INTRODUCTION 



An important problem arising in the study of supersymmetric models is the occurrence 
of non perturbative dynamical symmetry breaking |l| . The problem can be studied in 1 + 1 
dimensional lattice models where numerical tools are more effective and it should be easier 
to obtain a definite answer to the relevant questions. 

In particular, the simplest theoretical laboratory is in our opinion the N=l Wess-Zumino 
model that involves chiral superfields only and does not present complications related to 
gauge invariance. Rigorous results in the continuum case can be found in |2j. On the lattice, 
accurate numerical results are available pll^j, but nevertheless a clean determination of the 
supersymmetry breaking transition remains rather elusive. 

The numerical simulations in Ref. j^J were performed using the Green Function Monte 
Carlo (GFMC) algorithm and strong-coupling expansions. The physics of the model is fully 
determined by a single function of the scalar field, called in the following prepotential V(ip). 
All the results for the model with cubic prepotential indicated unbroken supersymmetry. 
Dynamical supersymmetry breaking in the model with quadratic prepotential V = A2<p 2 + A 
was studied performing numerical simulations along a line of constant A2, and confirmed the 
existence of two phases: a phase of broken SUSY with unbroken discrete Z2 at high Ao and a 
phase of unbroken SUSY with broken Z 2 at low Ao, separated by a single phase transition. We 
also studied the approach to the continuum limit in the model with quadratic prepotential 
performing numerical simulations along a 1-loop renormalization group trajectory in the 
phase of broken supersymmetry. 

The aim of the present paper is to propose a new approach to the study of the transition 
point. The method is based on the calculation of rigorous lower bounds on the ground state 
energy density in the infinite-lattice limit. Such bounds are useful in the discussion of the 
supersymmetry phase diagram as follows from a rather simple argument that we sketch here 
and expand later in the paper. The lattice version of the Wess-Zumino model conserves 
enough supersymmetry to prove that the ground state has a non negative energy density 
p > 0, as its continuum limit. Moreover the ground state is supersymmetric if and only 
if p = 0, whereas it breaks (dynamically) supersymmetry if p > 0. Therefore, if an exact 
positive lower bound plb is found with < plb < P, we can claim that supersymmetry is 
broken. 
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In the paper we find a sequence p^ of exact lower bounds representing the ground state 
energy densities of modified lattice Hamiltonians describing a cluster of L sites, which can be 
computed numerically. We stress that p^ is an exact bound for each L and moreover that 
it is consistent in the sense that - p from below for L — > oo. This features will permit 
to draw conclusions on the critical values of the coupling constants in the infinite-volume 
limit. 

The plan of the paper is the following: in Sec. Owe review the lattice formulation of the 
N=l Wess-Zumino model. In Sec. II I II we derive the lower bound. In Sec. IIVI we review the 
available results on the supersymmetric transition of the model and discuss our new results. 
Conclusions are summarized in Sec. Finally, in Appendix [X] we discuss the bound in 
the free limit and in Appendix El we collect a few useful technical details on the simulation 
algorithm used for the actual evaluation of the bound. 



II. THE N=l WESS-ZUMINO MODEL ON THE LATTICE 



The most general SUSY algebra in two dimensions has N left-handed fermionic gener- 
ators {Ql}a=i,...,n and N right-handed fermionic generators {Qr}a=i...n anc ^ * s denoted 
by (A, N). The bosonic generators are the components of the two-momentum (P°, P 1 ) and 
central charges T AB . In the (A, A) = (1,1) case we denote Qi j2 = Q\ ± Q\ and find 
{Q a , Qb] = 2{H 1 + Pa 1 + Ta 3 ) ab , where a* are the Pauli matrices, (P°, P 1 ) = (H, P) and 
T = T . The Wess-Zumino model realizes the above algebra on a real chiral multiplet with 
a real scalar component (p and a Majorana fermion with components ipip. The supercharges 
are 



Qi,2 = dx 



pry, - ( ||±1%>) ) 



(2.1) 



where p(x) is the momentum operator conjugate to tp(x) and V((p(x)) is an arbitrary function 
called prepotential in the following. 

On the lattice there are no continuous translation and we can only preserve a SUSY 
subalgebra [5|. We pick one of the supercharges, say Q x , build a discretized version Q L and 
finally define the lattice Hamiltonian to be H = Q\. 

The explicit lattice model is built by considering a spatial lattice with L sites and open 
boundary conditions. On each site we place a real scalar field (p n together with its conjugate 
momentum p n such that [p n ,y?m] = —iS n ,m- The associated fermion is a Majorana fermion 
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ip a , n with a = 1,2 and {ip a , n , ipb,m} = S ayb 5 n . m , i>\ n = ip a , n . The discretized supercharge is 



n=l 



(2.2) 



Following we replace the two Majorana fermion operators with a single Dirac operator c 
satisfying canonical anticommutation rules, i.e., {c n ,c m } = 0, {c^c^} = S ntm : 

V>{l,2},n = ^ 2 ; .„ (Ci±ic n ). (2.3) 

The Hamiltonian is H = Hb + -£^f with 

«» = t l^ + K ^T'" 1 + ' / ^>) 2 }- < 2 - 4) 

«f = £ |-i(4c„ +1 + fc.c.) + (-l)"V( Vn ) (c*c„ - i] I . 

n=l \ / ^ 

It conserves the total fermion number 

L 

N f = Y,N n , N n = cic n , (2.5) 

n=l 

and can be examined in each sector with definite Nf separately. In this work, we shall always 
consider the half-filled sector where the ground state is expected to lie. We remark that in 
this sector there exists a particle-hole symmetry which will play a role in the following. 

The relevant quantity for our analysis is the ground state energy density p evaluated on 
the infinite lattice 

p=}jm-°U.. (2.6) 

L^OO Li 

It can be used to tell between the two phases of the model: supersymmetric with p = 
or broken with p > 0. In the next Section we shall obtain rigorous lower bounds on p and 
exploit them to determine the phase. 

III. A FAMILY OF LOWER BOUNDS ON THE ENERGY DENSITY 
A. Derivation of the bounds 



Given a translation-invariant Hamiltonian H on the lattice Z it is possible to obtain lower 
bounds on its ground state energy density from a cluster decomposition of H. Briefly, given 

4 



a suitable finite sublattice A C Z d , it is possible to introduce a modified Hamiltonian H 
restricted to A such that its energy density p A bounds p from below. The difference between 
H and H is almost trivial and amounts to a simple rescaling of its coupling constants. This 
simple fact was already noted in 1951 by Anderson []| and has been recently improved and 
exploited in purely fermionic extended Hubbard models The derivation is quite general 
and applies to generic models with finite range interactions as we now review. 

As the sublattice A invades Z, the bound p A is expected to improve its accuracy with 
Pa — > p in the infinite A limit. This property, that we call consistency in the following, will 
play a crucial role in the analysis of the Wess-Zumino model and, for this reason, we shall 
discuss it in some details. 

As a preliminary discussion, it is expedient to make a general remark. Let H be any 
translation-invariant Hamiltonian on the d- dimensional lattice Z d , the only restriction on H 
being that of a finite-range interaction between non overlapping subsets. For instance, H 
could include nearest-neighbor interaction terms, next-to-nearest-neighbor couplings, and so 
on until some finite distance. 

Let A be any finite subset of the lattice, H A the restriction of H to A and Eq(A) the 
energy of its ground state |^o,a)- We define 

^def^^A)^ (31) 



and show that 



where for each local observable O, 



A->oo |A| 



(O) = lim (Vq,a|Wo,a) = lim (0) A . (3.3) 

A— >oo A— >oo 

The expectation value (O) means the average of O in the infinite- volume ground state of 
H . The proof of Eq. ()3.2|) is easy. Since (H\) = limA'(ifA)A' > E (A), we have the trivial 
inequality 

lim rrT (#a) > lim ^jyp = p. (3.4) 

A^oo | TV I A^oo | TV | 

We'll show now that the opposite inequality holds, and consequently that Eq. (J3.2)) is true. 

Due to the finite range of the interaction between non overlapping subsets of Z d , when 
A' is a sufficiently large region containing the bounded subset A in its interior, we have 

if A' = Ha + #a'\a + H dA , (3.5) 



where <9A is the boundary of A. Consequently, 

E (A') = E (A) + E (A'\A) + 0(dA). (3.6) 

Now, 

E Q (A') = (H A ,) A , = (H A + H A r\ A + Hq A ) A i > (H A ) A . + E (A'\A) + (H dA ) A ,, (3.7) 
and therefore, by exploiting Eq. ()3.6|) . we get the inequality 

(H A ) < E (A) + O(dA), (3.8) 

when A' — ► oo. Therefore 

lim ±(H A )< hm ^ = p. (3.9) 

A^oo A A^oo |A| 

Coming back to the lower bound or the ground state energy density, in order to illustrate the 
method we postpone the discussion of the Wess-Zumino model and consider first a generic 
lattice Hamiltonian with only on-site and nearest neighbor interactions 

oo 

H= £ (H^+hW), (3.10) 



i.e., Hn has an explicit form that does not depend on the base site n. We denote by H. 



where involves operators at site n, n + 1, . . . , n + k. We assume translation invariance, 

s,N 

the restriction of if to a (connected) sublattice starting at site s and ending at s + N — 1, 
i.e., a cluster of sites starting at s. We define it precisely as 

s+N-l s+N-2 
n=s n=s 

To derive the bound it is convenient to define the following cluster Hamiltonian where the 
relative weight of the on site and nearest neighbor terms has been modified 

s+L-l T s+L-2 

Hs,L= E + E H n } - ( 3 - 12 ) 

n=s n=s 

We now superpose several copies of H s l in order to reconstruct the initial Hamiltonian. We 
begin with a fixed lattice of A^ sites and define accordingly 

N-L+l 

H(n,l )= J- H n>L . (3.13) 

n=l 



We can expand the sum and rearrange terms obtaining 



(N-L+l N-L \ 

£ H^ + ^H^A+B, (3.14) 
n=L n=L J 

= LH L 

,N-2L+2 

where B is a boundary term spanning a number of sites growing like L, but independent 
on N, located at the left and right ends of the lattice. The above results is made intuitive 
in Fig. E] where we show the obvious origin of the factor L/(L — 1) in the nearest-neighbor 
terms of H u> l- 

It is now convenient to introduce 

Eq N ' L ^ = ground state energy of H^ N ' L \ (3.15) 
= ground state energy of H, ^, (3.16) 
p (L) = e^/L. (3.17) 

where in Eq. (13.16(1 we have not specified the initial site due to translation invariance. From 
Eq. 1)3.13(1 and convexity we deduce the inequality 

Ef' L) >(N -L + 1)E^ L \ (3.18) 

and dividing by LN 

Ell > ^£±V> (3.19) 

Let p* be the limit of the l.h.s. as iV — > oo with fixed L. If we prove that p > p* we obtain 
the desired lower bound 

p > p {L) . (3.20) 
To show this let us rewrite Eq. (J3.14j) in the form 

Hl,N-2L+2 H( n ' l ) B 



(3.21) 



N LN LN 

Taking the expectation value (•) with respect to the ground state of the original Hamiltonian 
on the infinite lattice we find 

(H L ,N-2L+2) > E^ N,L) _ (B)^ 

N ~ LN LN' 1 ' ' 

In the limit N — > oo, the contribution of the boundary terms vanishes because it does not 
depend on N due to the translation invariance of (■). The l.h.s. of Eq. ()3.22() is p in the 



limit, by the general remark we made at the beginning and we thus find p > p* completing 
the proof of Eq. (J3~2U|) since p > p* > p (L) . 

The bound is also consistent in the sense that 

lim p (i) = p, (3.23) 

L— >oo 

as it can be shown by an argument similar to the above one. To this aim, we write 

n=s 

and we take now the expectation value with respect to the vacuum of H s ,l m the infinite- 
lattice limit L — > oo. Since 

lim p^ = lim T (H 0>L ), (3.25) 

L— >oo L^oo J_, 

by the general remark we made at the beginning, we obtain from Eq. (|3.24j) the inequality 

lim p {L) > p, (3.26) 

L^oo 

that combined with p^ < p gives Eq. ()3.23|) . 

In the case of the Wess-Zumino model it is necessary to include also next-to-nearest 
couplings, but the procedure is identical and we just sketch its main features. We write 

oo 

H= E^f + ^+tff), (3-27) 

n=— oo 

(k) 

where Hn involves operators at site n, n + 1, . . . , n + k. We have explicitly 

#f = \vl + \pl + \v{<p n ? + (-l) n V'(<p n ) Uc n - , (3.28) 
= \v{ip n )ip n+1 - hp n V{ip n+1 ) - l -{clc n+1 + h.c), (3.29) 

#f = -\^n+2. (3.30) 



Again, we define 

H s ,l =E^ + I7lE ^ + J^- 2 £ R n ] - (3-31] 



-L-l T s+L-2 T s+L-3 



Due to the presence of the term (—1)™ in the fermion-boson coupling in Eq. (J3.28|) the 
Hamiltonians H s ,l for even and odd s are different, but they are related by the particle- hole 
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symmetry, as we remarked before. Therefore they have the same ground state energy and 
the previous discussion applies. Still, there is a translation invariance on our lattice, but 
only under shifts by two lattice spacings. In our staggered fermion formulation, it is the 
remnant of the continuum translation invariance. As before, we prove for 

= ground state energy of —H, L , (3.32) 

hi 

the relations 

P>P {L \ with limp (i) =p. (3.33) 
B. Relevance to the problem of supersymmetry breaking 

We now explain the way we exploit the sequence of bounds p^ (indexed by the cluster 
size L) to determine the phase at a particular point in the coupling constant space. To this 
aim, we compute numerically p^ L ' at various values of the cluster size L. If we find > 
for some L, we can immediately conclude that we are in the broken phase. If, on the other 
hand, we find a negative bound we cannot conclude in which phase we are. However, we 
know that p^ — ► p for L — > oo and the study of p^ as a function of both L and the coupling 
constants permit the identification of the phase in all cases. We shall discuss this strategy 
in more details and for the particular case of the Wess-Zumino model in the next Section 
dedicated to the presentation of our numerical results. 

We remark that the calculation of p^ L > is numerically feasible because it requires to 
determine the ground state energy of a Hamiltonian quite similar to H and defined on 
a finite lattice with L sites. In particular, it is not necessary to deal with the difficult 
problem of computing any infinite-size limit. In other words, we are able to transfer the 
information obtained in finite volume to a result holding on the infinite lattice. As a simple 
semi-analytical illustration of the technique, we report in Appendix El the evaluation of the 
bound in the free model with zero prepotential. 

IV. NUMERICAL RESULTS 

To test the effectiveness of the proposed bound and its relevance to the problem of locating 
the supersymmetric transition in the Wess-Zumino model, we study in details the case of a 
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quadratic prepotential 

V(<p) = W + Ac, (4.1) 

and discuss the dependence of p on A at fixed A 2 . Indeed, in the continuum, for a given 
A 2 , a general argument by Witten |9| and rigorous results in constructive quantum field 
theory [3] suggest the existence of a negative number Aq such that p(Ao) is positive when 
Ao > Aq and it vanishes for Ao < Aq. Such number Aq is the value of Ao at which dynamical 
supersymmetry breaking occurs. We draw in Fig. |2] a reasonable qualitative pattern of 
the curves representing p( L )(Ao). We see that a single zero is expected in p^(Ao) at some 
Ao — Ao(£). Since ]imj J -^ 00 p^ L ' = p, we expect that A (£) — > Aq for L — > oo allowing for a 
determination of the critical coupling Aq. 

We stress that the continuum limit of the model is obtained by following a Renormal- 
ization Group trajectory that, in particular, requires the limit A 2 — > 0. This step has been 
partWly — shed in Rrf . flashed in the next paragraph, but will not be pur- 
sued in the present paper. Here, we consider the lattice Wess-Zumino model at intermediate 
couplings and, we stress again, in the infinite-lattice limit. 



A. Review of GFMC results 

Here we give a brief review of our results in Ref. jjj, in order to compare them with 
the results in the present paper. In Ref. the analysis of supersymmetry breaking in a 
class of two dimensional lattice Wess-Zumino model (with open boundary conditions) was 
performed using the GFMC algorithm. This method computes a numerical representation of 
the ground state on a finite lattice with L sites in terms of the states carried by an ensemble 
of K walkers. The need to extrapolate to infinite L and K is the main source of systematic 
error. As an example of odd prepotential the case V = <p 3 was investigated, measuring 
the ground-state energy and the supersymmetric Ward identity. Both of them gave a very 
convincing evidence for unbroken SUSY. In Ref. (^J, in order to study the supersymmetry 
breaking pattern, the more interesting case of even prepotential (|4.1|) was investigated. In 
this case the model enjoys an approximate Z 2 symmetry which corresponds to the symmetry 
under the transformations <p n — > — <p n , Xn Xn- The theoretical expectation is that, for a 
fixed value of A 2 , at high Ao the model is in a phase of broken SUSY and unbroken Z 2 ; at 
low Aq it is in a phase of unbroken SUSY and broken Z 2 . The case A 2 = 0.5 was investigated 
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in detail. 

The usual technique for the study of a phase transition is the crossing method applied to 
the Binder cumulant, B, with a sensible choice of magnetization (which is not completely 
trivial, since our model is neither ferromagnetic nor antiferromagnetic and it doesn't enjoy 
translation symmetry). The crossing method consists in plotting B vs. Ao for several values 
of L. The crossing point Ag r (Z/i, L2), determined by the condition 

B(\Z,L 1 ) = B(\%,L 2 ) 

is an estimator of Aq. The convergence is dominated by the critical exponent u of the 
correlation length and by the critical exponent uj of the leading corrections to scaling 

A c r (L 1; L 2 ) = \ ( i ) + 0{L-r 1/ \Lr- 1,v ); 

we expect the phase transition we are studying to be in the Ising universality class, for which 
v — 1 and lu = 2, and therefore we expect fast convergence \™ — > Aq. The results [3| indicate 
Ao,gfmc = -0-48 ±0.01. 

It is possible to study the phase transition by looking at the connected correlation function 
Gd = (Lp n ip m ) c , averaged over all n,m pairs with \m — n\ = d, excluding pairs for which m 
or n is close to the border. Even and odd d may correspond to different physical channels. 
In Ref. (2! Gd is fitted to the form exp[— a\ — a 2 d + a 3 /(d + 10)], separately for even and 
odd d. In the broken phase, we have small but nonzero a 2 , and we observe equivalence of 
the even- and odd-<i channels, while in the unbroken phase, a 2 is larger, and the even- and 
odd-d channels are different The difference between the two phases is apparent, e.g., in 
the plot of a 2 vs. A . The data presented in Ref. 0] confirms the quoted value of Ag 1 . 



B. Results for the lower bound p^(Ao) 

We consider the quadratic prepotential Eq. (|4.1j) at the fixed value A2 = 0.5. As discussed 
above, the properties of the bound p^(Ao) guarantee that for large enough L it must have 
a single zero Aq(L) converging to Ao as L — > 00. In any case for each L we can claim that 

1 Other methods to estimate Aq were discussed in Ref. but with quite larger errors. In particular a 
study of the bosonic fields effective potential suggested Aq ~ —0.40, whereas an extrapolation of Eq/ L to 
infinite K and L give Aq ~ —0.53, with a rather large uncertainty. 
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Ag > Xq(L). To obtain a numerical estimate of p( L >(\o) we exploit the so-called worldline 
path integral (WLPI) algorithm. We discuss it in some details in Appendix |B] and we just 
recall here its basic features in order to introduce the systematic errors of the analysis. The 
WLPI algorithm computes numerically the quantity 

^.D^ wfi^a , (4 . 2) 

L Tr{(e T^i e T^)T} 

where the Hamiltonian for a cluster of L sites (I3.31J1 is written H — Hi + H 2 by separating 
in a convenient way the various bosonic and fermionic operators in the subhamiltonians Hi 
and if 2- The desired lower bound is obtained by the double extrapolation 

p {L) = lim lim p {L) {p,T), (4.3) 

with polynomial convergence ~ 1/T in T and exponential in (3. Numerically, we determined 
p( L \(3,T) for various values of (3 and T and a set of Ao that should include the transition 
point, at least according to the GFMC results. 

Our basic numerical results are shown in Figs. 0HZ1 where we plot the function p^(P,T) 
for L = 6-18, various (3 and T = 50, 100, 150. In order to avoid fermionic sign problems we 
need L = 4k + 2. In this case, at half filling, when a fermion hops through the lattice it passes 
over L/2 — 1 = 2k fermions, i.e., an even number giving no signs in the partition function. 
The extrapolation to T — * oo is quadratic in 1/T and a sample case is shown in Fig. |HJ The 
results are shown in Fig. El The chosen values of (3 are such that the curves corresponding 
to the highest (3 values can be taken as representatives of the (3 — > oo limit. This is rather 
safe due to the exponential convergence in (3. For all cluster sizes, we see that the energy 
lower bound behaves as expected: it is positive around Ao = and decreases as Ao moves 
to the left. At a certain unique point Ag(L), the bound vanishes and remains negative for 
A < Aq(L). This means that supersymmetry breaking can be excluded for A > min^ Aq(£). 
Also, consistency of the bound means that Ag(L) must converge to the infinite- volume critical 
point as L — > oo. Since the difference between the exact Hamiltonian and the one used to 
derive the bound is 0(1/L), we can fit \q{L) with a polynomial in 1/L. This is shown in 
Fig. El where we also show the GFMC result. The best fit with a parabolic function gives 
Aq = —0.49 ± 0.06 quite in agreement with the previous Aq GFMC = —0.48 ± 0.01. 

The somewhat large error could of course be reduced by improving the (3 and T extrap- 
olations and with additional statistics; however, we do not pursue further the numerical 
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analysis, since the main aim of the present paper has been to show the validity of the new 
proposal for the identification of the infinite-lattice critical point, by a rather non-standard 
method. 

As a concluding remark, we briefly discuss the cubic prepotential V(ip) = |y? 3 , where 
general arguments predict no supersymmetry breaking; this has been confirmed by the 
GFMC analysis jsj]. In Fig. ^2 we show the lower bound computed on the clusters with 
L = 6, 10, 14 and 18 sites with (3 = 10, 12, 14 and three values T = 50, 100, 150. We also 
show the results of a quadratic extrapolation in XjT at fixed (3. At a rather qualitative level 
we remark two features of the Figure. First, the bound is definitely negative for all L; this 
is consistent with our expectations because the infinite-volume energy density vanishes in 
this case. Second, the values at large (3 approach zero as L increases as follows from the 
discussed consistency properties of the bound. 



V. CONCLUSIONS 



In the present paper we have studied the phase diagram of the lattice N=l Wess-Zumino 
model with a particular scalar potential for which we expect dynamical supersymmetry 
breaking to occur. As it is known, the model can be put on the lattice with a certain 
amount of supersymmetry left unbroken by the space-time discretization. As an important 
consequence, the lattice model shares with its continuum counterpart the important property 
of having a non negative energy density p > 0. Also the vanishing of p is a necessary and 
sufficient signal for unbroken supersymmetry, whereas p > would imply a dynamically 
broken phase, that is expected to exist on the infinite lattice. All these facts trigger the 
interest toward a tight lower bound on p holding, we repeat, on the infinite lattice. We have 
proposed precisely a family of such bounds p( L > that can be computed in terms of the energy 
density of a suitable (non supersymmetric) local model defined on a cluster of L sites. The 
family is consistent, i.e., p^ — > p as L — > oo. 

The analysis of p( L ' as function of the cluster size L and the coupling constants permits 
to derive bounds on the critical values of the couplings where dynamical supersymmetry 
breaking occurs. Also, extrapolation in the cluster size allows in principle to locate the 
transition in a rather straightforward way. 

To confirm the feasibility of the procedure we have computed numerically p^ L \ with good 
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agreement with existing previous studies based on extensive simulations with the Green 
Function Monte Carlo method. We thus conclude that the proposed approach is useful as 
an independent method to locate the supersymmetric transition. 

As a final comment, we would like to emphasize what seems to us a general interesting 
point. The condition p > in (globally) supersymmetric models is obviously and directly 
related to the issue of symmetry breaking. Usually, this simple fact is not exploited in non- 
perturbative approaches based on standard simulations in the Lagrangian framework where 
the focus is mainly on correlation functions to check for instance Ward identities. This is 
definitely not true in the (lattice) Hamiltonian formalism where the ground state energy is 
typically the simplest observable that can be considered, at least in finite volume, and the 
vanishing of p is the most natural monitor of supersymmetry. As we have shown, analytical 
information on p, like rigorous bounds, can be immediately translated into a valid tool to 
make the numerical investigation most effective. 
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APPENDIX A: EVALUATION OF THE BOUND IN THE FREE MODEL 



The Wess-Zumino model with linear prepotential V(ip) = pip is free with bosons decou- 
pled from fermions. The Hamiltonian H^l can be written 

l „ L L 



= 5 Eft 2 + 5 E <W, + E <M<* 



where the matrices and have elements 



fti=i 



-Jzb \i-3\ = 2 
\ + p 2 i = j 
else 



M {F) 



m |i-j| = i 



p(-lYt = j 

else 



(Al) 



(A2) 



They can be diagonalized and we denote by Ai < A 2 < • • • < the eigenvalues of and 
by Luf, ■ ■ ■ ,u\ those of M^ B \ The ground state factorizes in the tensor product of a bosonic 
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and a fermionic ground state. Their energies are given respectively by 

4=1 i=l 

where we assumed L to be even. These expressions are simply the ground state of a multi- 
dimensional quantum harmonic oscillator and that of a system of free (non-relativistic) 
fermions at half filling. The bound on the ground state energy density is 

P iL) = j(E { B) +E^). (A4) 

It can be evaluated numerically without difficulty for large L and its plot is shown in Fig. (J3J) 
for various fi. The model with linear prepotential is known to respect supersymmetry 
(p = 0). We thus expect that p^ < for all L and also that p^ — > as L — > oo 
with a leading correction ~ 1/L. These expectations are fully in agreement with the figure. 



APPENDIX B: DESCRIPTION OF THE SIMULATION ALGORITHM 

The world-line path integral computes the expectation value of an operator A at inverse 
temperature f3 as Tr(A e _/3H )/Tr e~^ H . As usual the trace is computed by splitting e~^ H 
and inserting with time spacing e complete sets of intermediate states that are sampled 
stochastically. 

In more details, we perform a first splitting separating the bosonic terms and the fermionic 
ones (we omit the base site because of spatial translation invariance) 

h& = hf + (Bl) 



where 



h ( B L) = £ (k + \n<P»? + lA) + (B2) 

n=0 ^ ' 

1 L L ~ 2 

2 J^TJ 22{V(<p n )<Pn+l - <PnV(ip n +l)) + (B3) 



71=0 

L-3 



4 L-2 

71=0 

and the remaining fermionic terms define hp as 



£ (B4) 



h { pl = -~-^c^ x + h.c.) + +c»(-i) b v( V 7„) (iv„ - 
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+C„ + i(-l) n+ V(^ n+1 ) (N n+1 - , (B5) 

where ( n is 1 on sites and L — 1 and 1/2 elsewhere. 
The calculation of Tre _/3h is written as usual 

Tre- ph = lim {^ {0) \e- £hB e~ £hF \^ T - l) ) ■ ■ ■ {^ l) \e~ £hB e- £hF \^) (B6) 

^ V'(o),v>( 1 ),...,^( T - 1 ) 

where Te = (3 and l?/?^) denotes a complete set of states at the kth. time slice. We work in 
a basis of states where the bosonic field ip and the fermion occupation numbers are diagonal 
and write therefore 

\^ k) ) = \^ {k \n ik) ) (B7) 

where (p^ k \ rS k "> are the eigenvalues of the above operators (vectors of L numbers, one for 
each site). The trace can be written in the following form 

T— 1 

Tr e^ h = lim £ I J[ Vip {t) e - s ^ +l ° sW ^ n \ (B8) 

n(0),n( 1 ),...,n( r - 1 ) *=0 

where the action Sb(^p) is 

= 53 2^(v».*+i ~ v™.*) 2 + £ 53 v * 

/ j 1 1 L L_2 

V * = 53 ( 2 y ^ n '*) 2 + 4^™-* + 2 L^T 53( y (^^ n+1 ~ Vn^(^+i))+ 

n \ n=0 

1 L L ~ 3 \ 

-| 53 <^«+2 1 ( B9 ) 

n=0 / 

and the weight W is the product 

W(<p,n) = (n^le^ln^-V) ■ ■ ■ (nW\e- ehF \n®). (BIO) 

In order to build a definite sampling procedure, it is convenient to elaborate further the 
explicit expression of W. This can be achieved by remarking that the Hamiltonian hp can 
be further split as 

hf = E ft S + E ft S- < B11 > 

n even n odd 

The even terms commute among themselves as the odd ones do. It is then convenient to 
put fermions on a square lattice of plaquettes as shown in Fig. H21 where the rows denote 
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alternatively evolution with respect to the even/odd terms. The weight W is obtained 
by multiplying the matrix element of exp(— ehp) for each shaded plaquette with lower-left 
corner satisfying (x + t) mod 2 = 0. 

If we denote the initial state of such a plaquette by \N n , N n+ i) and define t — L/{L — 1), 
a = C n (— l) n V'((f n ), b = ( n+ i(—l) n+1 V (f n+ i) , we have the basic results 

e- £HF < n \0,0) = e^ {a+b) \0,0), e- £HF '"\l,l) = e~^ a+b) \0, 0), 

e - efl >.«|l 0) = (coshep- -sinhep)|l,0) + — sinheplO, 1), 

p 2p 

e -eja>,„| i\ = (coshep + -sinhep)|0, 1) + — sinhep|l,0), 

p 2p 



where 



1 



t = j— i , r = —i p = -V^+¥ (B12) 



The full weight W is obtained by multiplying the suitable coefficient on the r.h.s. (depending 
on the final state on the plaquette) over all plaquettes 2 . 

After this preliminary discussion of the explicit form of the lattice partition function, we 
turn to the sampling algorithm and also discuss how energy is measured. 



1. Monte Carlo sampling 

The configuration update is obtained by first updating the bosonic fields in the fixed 
fermion background and then updating fermions with fixed bosonic fields. The bosonic 
update is a multi-hit Metropolis sweep. The fermionic update is an accept-reject step per- 
formed on admissible local changes of the fermion world-lines. These are described by the 
transformation |0, 1) — > |0, 1) into |1,0) — > |1,0) and vice versa on all the plaquettes with 
(x + t)mod2 = 1. A very clean detailed discussion of such alg orithms in the context of 
models with phonon-electron coupling can be found in Ref. |lol ] (see also Ref. 0| for the 
specific case of the Wess-Zumino model). 

2 Notice that W > because t > 0, p > and also 

coshepi -sinh e/ 9= i ( 1 + - ) e ep + i ( 1 - - ) e~ £p > (B13) 
P 2 V PJ 2 V PJ 

due to \r/p\ < 1. 
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2. Energy measurement 



Given the set of sampled configurations, we now explain how to measure the average 
energy, i.e., the quantity Tr(H e _/3H )/Tr e~^ H that we determine, of course, at the discretized 
level. For notational convenience we shall denote it by (H)p = (hs)^ + {h F )p. 

The calculation of (hs)^ is trivial apart from the p 2 part. This term is evaluated by the 
virial theorem 

-2\ v,,.. i . — \ ( B14 ) 



lim {p n )p 

p— >oo 



lim 



if., 



The r.h.s. does not present any problem because the fermionic terms appearing in dh/d(p n 
are diagonal in the chosen occupation numbers basis. 

The evaluation of (hp)p can be reduced to the calculation, for a given configuration, of 
the following basic matrix elements [lol ] 



(f\e- £hp h F \i) {f\h F e- £hF \i) 



(f\e-^\i) 



(B15) 



where \i) and |/) can be |00), |11), 1 01) , 1 10) . For transitions without hopping these read 



-ehf 



-ehi 



\n, n) 

and for transitions with hopping 



(n, n 


h F e~ £h 


F 


n, n) 


(n, n 




n, n) 



in ■ 1>) ( n- i 



(B16) 



and 



(l,0|e- eA|? /i F |l,0) 


a 


-b 


t k 


(l,0|e- £ ^|l,0) 




2 


2k^ 


(0, l\e- £hF h F \l,0) 




a + b 


t k_ 


(0,l|e- £ ^|l,0) 




2 




(0,l|e" £ ^/i F |0,l) 




a + b 


t k 


(0,l|e- £ ^|0,l) 




2 


2k + 


(l,0|e- £ ^/i F |0,l) 


a 


-b 


t k + 


(l,0|e- eh *|0,l) 




2 


2 V 


(l,0|/i F e- £ ^|l,0) 


a 


-b 


t ko 


(l,0|e- £ ^|l,0) 




2 


2k^ 


(0,l|/i F e- £ ^|l,0) 


a 


-b 


t k + 


(0,l|e" £ ^|l,0) 




2 


2 V 


(0,l|/i F e~ £ ^|0,l) 




a + b 


t ko 


(0,l|e- £ ^|0,l) 




2 


2k + 


(l,0|/i F e- £ ^|0,l) 




a + b 


t k_ 


(l,0|e- £ ^|0,l) 




2 


~2~ko 
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where 

t T 

k = —sinep, k± = cosh ep ± - sinh ep (B17) 



[1] A. Feo, Nucl. Phys. [PS] 119, 198 (2003). 

[2] A. Jaffe and A. Lesniewski, Supersymmetric quantum fields and infinite dimensional analysis, 
lectures published in Cargese Summer Inst. (1987) and references therein. 

[3] M. Beccaria, M. Campostrini, A. Feo, Nucl. Phys. [PS] 106, 944 (2002), 119, 891 (2003); M. 
Beccaria, M. Campostrini, A. Feo, Nucl. Phys. [PS] 129-130, 874 (2004); M. Beccaria, M. 
Campostrini, A. Feo, hep-lat/0402007, to appear on Physical Review D. 

[4] M. Beccaria and C. Rampino, Phys. Rev. D67, 127701 (2003). 

[5] S. Elitzur, E. Rabinovici, A. Schwimmer, Phys. Lett. B119, 165 (1982); S. Elitzur and 

A. Schwimmer, Nucl. Phys. B 226, 109 (1983). 
[6] J. Ranft, A. Schiller, Phys. Lett. B138, 166 (1984); 
[7] P. W. Anderson, Phys. Rev. 83, 1260 (1951). 

[8] H. Nishimori and Y. Ozeki, J. Phys. Soc. Jpn. 58, 1027 (1989); R. Tarrach and R. Valent, 
Phys. Rev. B 41, 9611 (1990); H. Wischmann and E. Muller-Hartmann, Phys. Rev. B 43, 
8668 (1991); T. Wittmann and J. Stolze, Phys. Rev. B 48, 3479 (1993); J. de Boer and 
A. Schadschneider, Phys. Rev. Lett. 75, 4298 (1995); Z. Szabo, Phys. Rev. B59, 10007 (1999). 

[9] E. Witten, Nucl. Phys. B188, 513 (1981); E. Witten, Nucl. Phys. B202, 253 (1982). 
[10] J. E. Hirsch, R. L. Sugar, D. J. Scalapino, R. Blanckenbecler, Phys. Rev. B26, 5033 (1982). 



19 



20 



FIG. 1: Explanation of the factor L/(L — 1) in the nearest-neighbor terms of the modified cluster 
Hamiltonian Eq. (|3.12|) . The black dots are the sites of a cluster with 5 sites. Circles and the 
arcs denote on-site and nearest-neighbor interactions. If the cluster Hamiltonian is shifted over all 
positions of an infinite lattice and all the contributions are summed, the dashed on-site interactions 
appear with multiplicity 5 and the dashed arc with multiplicity 4. In the general case, the ratio is 
precisely L/(L — 1). 
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